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We exploit a relation between the mean number A/'m of minima of random Gaussian surfaces 
and extreme eigenvalues of random matrices to understand the critical behaviour of A/'m in the 
simplest glass-like transition occuring in a toy model of a single particle in A'^-dimensional random 
environment, with A'^ ^ 1. Varying the control parameter /j, through the critical value fic we analyse 
in detail how Mm{fJ,) drops from being exponentially large in the glassy phase to Afm{fJ,) ~ 1 on the 
other side of the transition. We also extract a subleading behaviour of A/'m (m) in both glassy and 
simple phases. The width Sji/ of the critical region is found to scale as A^"^'^"^ and inside that 
region Mm (a*) converges to a limiting shape expressed in terms of the Tracy- Widom distribution. 
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Understanding statistical structure of stationary 
points (minima, maxima and saddles) of random land- 
scapes and fields of various types is a rich problem of 
intrinsic current interest in various areas of pure and ap- 
plied mathematics [J [5] . It also keeps attracting steady 
interest in theoretical physics community over more than 
fifty years (SHE], with recent applications to statistical 
physics [9l [10] , string theory [11] and cosmology [121 US] . 
For a landscape described by a sufficiently smooth ran- 
dom function 'H(x) of N real variables x = (xi, xa?) 
the problem of counting all stationary points amounts 
to finding solutions of the simultaneous stationarity con- 
ditions dkH — for all k = 1,...,N, with dk standing 
for the partial derivative . Finding the total number 

^ioP " Id of stationary points in any spatial 

domain D amounts to knowing the corresponding density 
/9(x) essentially given by the Jacobian associated with 
the Hessian ^ T-L oi the random surface at stationary 
points. In particular, the mean value of such density is 
given by the so-called Kac-Rice formula 

P™W = {\^^t{dlMMY[5{dkH)^ (1) 

where 5{x) stands for the Dirac's 5— function and brack- 
ets here and henceforth denote the ensemble average. 
For a general random surface the problem of evaluat- 
ing the averages involving the modulus of the Jacobian 
is rather difficult and no efficient technique seems to be 
known to perform the task [M]. However, as was first no- 
ticed in [T' such calculation can be indeed completed for 
Gaussian fields 'H(x) such that their covariance structure 
{'H{xi)H{x2)) depends only on the Euclidean distance 
|xi — X2I and is therefore invariant under rotations. For 
such fields it turns out to be possible to reduce ([T]) to eval- 
uating the mean density of eigenvalues of the Gaussian 
Orthogonal Ensemble (GOE) of real N x N random ma- 



trices for which closed- form expression is well-known [15j . 
This observation and its further ramifications proved to 
be very useful in estimating the probability density of the 
highest maximum of such surfaces (2j , for obtaining de- 
tailed information about large- iV asymptotics of A^^^'' in 
spherical spin-glass models [TB] and in counting station- 
ary points of random superposition of eigenfunctions of 
Laplacian operator on high-dimensional manifolds [17 . 

Among various possible types of random Gaussian 
landscapes we will concentrate on arguably the simplest, 
yet nontrivial model 

N 

fc=l 

where /i > is the main control parameter and V{x) is a 
random mean-zero Gaussian-distributed field character- 
ized by a particular (translational invariant) covariance 
structure: 

(nx)ny))=A^/(2^(x-yf) (3) 

where f{x) is any smooth function which in this paper 
we consider suitably decaying at infinity. 

Looking at ([2]) as a certain random energy surface, 
one can associate with it the equilibrium statistical me- 
chanical model of a single particle equilibrated by ther- 
mal forces as described by the corresponding equilibrium 
Boltzmann-Gibbs measure p^(x) = e^^^'^^'' /Z{(3), with 
/3 = being the inverse temperature and Z{/3) being 
the associated partition function. The model can be then 
studied in the framework of the replica trick method, see 
[T5H^ , though nowadays a rigorous mathematical treat- 
ment of models of such type is possible as well, see e.g. 
|21| . The analysis reveals that for ji < fic — \J /"(O) 
the model exibits for A^ — > 00 a well-defined thermody- 
namic transition to a glassy phase with single-step broken 



2 



replica symmetry. In particular, in the latter phase con- 
fined to the lower-temperature region T < Tc{^) in the 
^jL~T plane the thermodynamic expectation value of the 
particle displacement / x^p^ (x) dx becomes T— inde- 
pendent, i.e. particle remains "frozen" below the tran- 
sition line. It is therefore conventional to consider such 
model as a toy system describing the glass-like freezing 
transition. The transition temperature Tc(/i) tends to 
zero when /i — !■ /ic, so that for /i > /i^ even for T — Q 
the system remains in the replica symmetric phase. On 
the other hand for ^ < /Xc the replica symmetry at zero 
temperature is broken. In this way the model provides 
a rather unique possibility of studying zero-temperature 
glass transition by varying the control parameter /i. 

As the statistical mechanics at T is expected to be 
dominated by minima in the energy landscape, the exis- 
tence of the zero-temperature phase transition suggests 
that the number of landscape minima Mm in the two 
phases should be qualitatively different. Namely, it is 
expected that Nm{^J') <x. exp{NYlrn{^j^)), with > 

for fi < fic when the random part of H dominates, but 
drops to a subexponential value for fi > fic when the 
harmonic part of H dominates and tends eventually to 
Afm — 1 for /i — > OO. In such a context, it is natural to 
call the quantity 'EmifJ') > the complexity of minima. 

Even for such a simple model ^ understanding the 
statistical properties of the number of minima AfmifJ-) 
is a difficult, and to large extent open problem. Some 
progress is possible if one concentrates on studying the 
simplest informative quantity, the mean value (AfmilJ-))- 
And even analysing that mean as a function of fi is not 
at all trivial, and is actually the subject of the present 
Letter. To that end let us recall that the original paper 
[7] concentrated on analysing even simpler quantity, the 
mean of the total number of stationary points Nf^^ (/i) , 
and on extracting the associated cumulative complexity 
Tit (m)- Further developments of the method proposed 
in [22], and applied to the present model in [53] allowed 
to evaluate the complexities associated with the mean 
number of stationary points with any extensive index / — 
aN ^ 1 with a > 0, where the index / is defined as 
the number of negative eigenvalues of the Hessian at the 
stationary point. One then can consider the formal limit 
a — )■ and in this way extract the supposed complexity 
of minima which turned out to be given by [23] 



in) 



2/i 



- when fi < fj,c (4) 



and TmifJ-) = for fi > fic- Note however that a Q 
can not distinguish genuine minima with / = from any 
other saddle-points with non-extensive < I N. 

The goal of this letter is to perform the accurate eval- 
uation of (Mm) directly from the first principles, and in 
particular to analyse in great detail the so-called criti- 
cal regime, i.e. the vicinity of the critical point fi = fic 



where the number of minima drops from its exponentially 
big value in the glassy phase fi < fic to a. subexponen- 
tial value for fj, > fic- Our starting point is the gen- 
eral Kac-Rice expression for the mean number of minima 
{^fm) = J p„i(x)d^x with p„(x) now given by (cfjl]) 



N 



p„(x) = ( |det (dl^n) I 9 (dl^n) n ^ (dk-H) 



(5) 



where the matrix Heaviside step-function 0{A) = 1 for 
positive definite matrices and zero otherwise, which en- 
sures that only strict minima are contributing to the 
counting. For the model in question ^ we can now 
follow the method of [7 and show that 



+ 00 



dt 



e-'^'-^KMizt) (6) 



where Zt = fi -\- fict and we introduced 

KNiz) = ( |det (z - Afo)| {z - Mo) ) ^^^^ (7) 

where the average now goes over the random ma- 
trix Mq taken from the so-called Gaussian Or- 
thogonal Ensemble (GOE) [T^ with the proba- 
bility density P(Mo) = Cn exp | - trM^ | . 

Using the 0{N) invariance of the GOE mea- 
sure, we can introduce the eigenvalues of Mq 
and in the standard way [15 1 find that K]\[{z) — 

z^^(2/i^/iV)^T KAr(y) with z and y related by z 

the normalization factor zn given by: 



ZN = Mc 



and where 



dXi... 



N N 

(2^)f llT{l+j/2) 
(8) 



N 



d\NY[\\^-\J\Y[iy-k 



i<j 



(9) 

In fact in the mathematical literature one frequently 
uses the "standardized" GOE defined as the ensemble 
of real symmetric N x N matrices M with the measure 
P(M) cx To that end define the partition 

function ZN{y) as 



Zn{v) 



(10) 



By definition, the probability that the maximal eigen- 
value of a standardized GOE matrix M is smaller than 
some value y is given by: 



^jv(Amax < y) 



Zn{v) 
Zn{<x}) 



(11) 
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Our key observation is that the function kjv (y) is related 



to the above cumulative distribution (11) as 



dZNjy) 
dy 



Ne 2 kN-iiy) 



(12) 



Moreover, a simple change of variables shows that 

N(N+1) 



ZNioo) = ZN 



. In this way we can ex- 
press the mean number of minima of our random energy 



surface as 



{Af„,) = (f) \lJ iMf^c) where 



iNifJ-ffJ-c 



+ 00 



(13) 



X— [PAr+i(A,„ax < v)] dy 

dy 



with PAr(Aniax < J/) defined in (11). It is appropriate to 



mention that a closely related formula appeared also in 
[l6] in the context of studying the extrema of the random 
energy surface for the so-called spherical model of spin 
glasses. 

Main utility of the above observation is in the fact 
that the behaviour of the cumulative distribution of the 
maximal eigenvalue of GOE matrices for large N was 
thoroughly studied in various regimes in recent years, 
starting from the famous work by Tracy and Widom [21] . 
In particular, most detailed large-deviation results were 
obtained in [25] for the left tail and [27] for the right tail 
for the regime N ^ 1 and y — s\/N with fixed s. For 
example, the left tail is given by 



d 
dy 



V(Amax < y)] 



for y < V2N (14) 



where y - V2N = 0{y/N) and 



N^tP^{s) - iV$i(s) + (j)i \nN + $2(s) 
108 ^ 3(3. 



s 

108 
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In 



+ Vs^ + 6 
3x/2 



(15) 



and explicit expressions for $i,2(s) and 01 are rather long 
and can be found in [25j . On the other hand, the right 
large deviation tail is given by [57] 

^ g-7VV+(s) 

-f- [Pw(A„,ax < y)] , 

dy 20F (-2 + s2)i/4 ^ yf=2+^ 



for s = 



/N 



and y > V2N where [25].[?7] 



1 - 



In 




(16) 



(17) 



These results allow us to get the behaviour of the mean 
number of minima by extracting the appropriate asymp- 



totic from the integral ( 13 1 for large N via the saddle 



point method. Using (16), we find for /i > /i^ the follow- 
ing equivalent: 



1 when II > fj.c 



(18) 



which is much more precise than just the statement of 
vanishing complexity. It turns out we just have on av- 
erage one single minimum not only for fj, ^ iic SiS one 
might nai vely expect, but immediately for all jj, > n^. 
Similarly, (14) implies that {Mm) ~ e'^™(^/'^'=) for fi < Hc 
where 



'TV 



(19) 
(20) 



with the leading term S (y) = — In (y) — ^ -|- 2y 
coinciding with the complexity of minima Eq. (|4| 
found in |23| and subleading terms given by Si (y) 
IV2{1 - yf'\ S2 (y) = -|23/4 (i _ yf^ and 



S3(y) 



-2 + 2y 
23 



1 



4 

+ gln([l-y]N'/^) + '-C{^l) 



y~ 



137 In 2 hiTT 
96 ^ ~Y~ 



(21) 
(22) 



It is easy to check that E (//) cx (/ic — A*)'^ when approach- 
ing the transition point as was already noticed in [23| . 
This should be contrasted with the cumulative complex- 
ity for all stationary points (instead of minima) vanishing 
as (/i) cx {lie — m)^ [7]- The difference is significant as 
it has implications for the width of the so-called "tran- 
sition region" |/i — nd where the two phases of the sys- 
tem become indistinguishable for large but finite N ^ 1. 
Here we argue that it is the complexity of minima which 
provides the correct scaling. Indeed, as shown in [20] the 
difference AF between the zero-temperature free energies 
F = — limT->.o Tin Z(/3) of the replica-symmetric solu- 
tion and one with broken replica symmetry is of the order 
of N{nc — 1^)^ close to the transition point. Thus ther- 
modynamically the two phases cease to be distinguish- 
able precisely in the same region when the leading term 

TV E in the lnA/'„i(^) becomes of the order of unity. 
One may consider this as an independent confirmation of 
the thermodynamic relevance of minima, rather than of 
the totality of stationary points for the glass transition. 
Moreover, we further see that for |/i — ^ A^^^/^ 

not only the leading term N E but all the sub- 

leading terms in ( 19 ) become simultaneously of order 
of unity. This fact implies that in the transition re- 
gion where (^j^ — 1^ iV^/^ = (5 is of the order of unity 
the mean number of minima tends to a N— independent 
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function of the scaling variable 6. To find that limiting 
shape we use the celebrated Tracy- Widom law ^24j for the 
probability of the maximal eigenvalue of the standardized 
GOE matrix: 




where J^i{x) is a special solution of the Painleve II equa- 



tion. Using this fact we can rewrite (|13|) as 
hN{x) 



(24) 

lnJ'i'(a;)(25) 



and = 2{N + I) + x ^^ +^ " . For large N, the inte- 
gral computed by the saddle point method is dominated 
by the vicinity of x* such that 




d 

dx 



\nF^{x*) 



(26) 



We thus get {Mm) ^ A/'((5) in the transition region = 
1 + SN^^ where the limiting shape J\f{S) is given by 



\ 



27r 



2J"i(a;*) e 



:*i5- 



(27) 



that the scaling function (27) which is one of the main 



results of this work matches smoothly the two regimes 



/i < cf ( 19 1 and /i > cf ( 18 1 as shown in the figure: 



1.10 

^J■/^|■c 



FIG. 1: Plot of the shape of the (log of the) mean number 
of minima {^^„^) for iV = 10000 (a) for /i < fic- {Mm) ~ e^" 
(blue dashed line) (b) in the transition region = 1+SN~^^^ 
(red solid line) (c) for fi > {Mm) ~ 1 (black solid line) 



In conclusion, we investigated in a great detail the 
mean number of minima of the random energy function 



in a toy model of the glass transition, and determined 
the precise scaling form of that number at the transi- 
tion region between the two phases. It turned out to 
be related to the Tracy- Widom distribution well-known 
in the random matrix theory. Note that such distribu- 
tion was also shown to describe fluctuations of transi- 
tion temperature in mean- field spin glasses [53]. We hope 
that similar methods could prove useful in the problem 
of understanding the ground states fluctuations of glassy 
systems, which attracts growing interest in recent years 
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